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Abstract 

A new class of bivariate distributions is introduced that extends the 
Generalized Marshall-Olkin distributions of Li and Pellerey (2011). Their 
dependence structure is studied through the analysis of the copula func¬ 
tions that they induce. These copulas, that include as special cases the 
Generalized Marshall-Olkin copulas and the Scale Mixture of Marshall- 
Olkin copulas (see Li, 2009), are obtained through suitable distortions 
of bivariate Archimedean copulas: this induces asymmetry, and the cor¬ 
responding Kendall’s tau as well as the tail dependence parameters are 
studied. 
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1 Introduction 

In their seminal paper (see Marshall and Olkin, 1967), the authors introduce 
the Marshall-Olkin distribution whose survival version is 


F{x\,X2) = exp{—Ai^i — X2X2 — A3 max(a;i,a;2)} 


xi,X2 > 0 , Ai,A2, A3 > 0. This is the distribution of a bivariate random vector 
{Ml, M2) with Ml = min(Ai,A3) and M2 = min(A2,A3), where Xi, X2, A3 
are three independent and exponentially distributed random variables. 

For an interpretation of the Marshall-Olkin distribution, consider a system with 
two components Ci and C2 (electronic elements, engines, credit obligors, life- 
insured married couples, etc.). The random variables Xi and X2 represent the 
arrival times of a shock causing the failure or default or death (depending on the 
case) of Cl and C2, respectively (idiosyncratic shocks), while A3 is the arrival 
time of a shock causing the simultaneous end of the lifetimes of both Ci and C2 
(systemick shock). 
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In order to generalize this model, the case in which the underlying random 
variables Xi, X2, X3 are assumed to have marginal distributions different from 
the exponential one has been considered in the existing literature. The most 
general results, in this direction, are those obtained in Li and Pellerey ( 2011 ), 
where no restriction is made on the marginal distributions of the random vari¬ 
ables Xi,X2, X3 while joint independence is again assumed. More precisely, the 
authors assume for (Xi, X2, X3) a joint survival distribution of type 

F{x,y,z) = exp(-(iLi(a;) + H2iy) + Hsiz))) ( 1 ) 

where the Hi are the cumulative hazard functions of the random variables Xi. 
The corresponding joint survival distribution of the random variables Mi and 
M2 is 

F{xi,X2) = exp{-Mi(xi) - H2{x 2) - H3{ma,x{xi,X2))} 

that is called Generalized Marshall-Olkin distribution. These distributions in¬ 
corporate, as special cases, the Marshall-Olkin type distributions introduced 
in Muliere and Scarsini ( 1987 ), the bivariate Weibull distributions introduced 
in Lu ( 1989 ) and the bivariate Pareto distributions introduced in Asimit et al. 
( 2010 ). The authors analyze the dependence structure implied by these distribu¬ 
tions introducing the corresponding copula functions, that they call Generalized 
Marshall-Olkin copulas. 

However, because of the independence among the original shocks (systemic 
and idiosyncratic), the dependence between Ci and C2 lifetimes is only given 
by the occurrence of the systemic shock. None the less, one could imagine 
concrete situations in which there is some additional dependence between Ci 
and C2 lifetimes given, for example, by some unobserved factors affecting all 
the original shocks (idiosyncratic and systemic) arrival times. 

Consider for example the case of a hrm interested in protecting itself against 
the failure of a production chain composed by a primary machine (or electronic 
component) in serial connection with two secondary machines in such a way 
that the failure of the primary machine causes the failure of all the productive 
system, while the failure of a secondary machine causes only the failure of the 
corresponding production line. Clearly the disease of all the three machines 
(representing the systemic and the idiosyncratic components of the system) is 
influenced by factors such as maintenance, good electric supply, and so on, and 
this fact induces dependence among them so that the disease of a production 
line can induce a change in the probability of failure of the primary machine. 
An insurance policy written in order to protect against lack in the production 
caused by the failure of all, or part, of the system, has to take into account the 
probability of failure of each machine (systemic and idiosyncratic effects) and 
the dependence among them. 

The above is a lucky case, in which the systemic and the idiosyncratic com¬ 
ponents are well identified and observable. This is not in general the case. 
Consider, for example, a life insurance contract written by a married couple 
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with a final payment or goods transfer made when both spouses have died: 
such contracts are thought in order to transfer money or goods to heirs or to 
guarantee some annuity until death, using, for example, the house property as 
collateral security (this is the case of reverse mortgages contracts). In such 
cases the contract expires when both spouses have died and their death can 
occur separately or simultaneously. Their simultaneous death can be caused by 
the occurrence of some systemic event (for example some “catastrophic” event 
like a car accident, an earthquake and so on) but also by something more re¬ 
lated to some dependence among their lifetimes idiosyncratic components: it is 
in fact well known that in very old people the death of one of the two spouses 
can induce in a very short period (so that the two events can be considered 
simultaneous) the death of the other. However, notice that, some catastrophic 
events causing the simultaneous death of the married couple, cannot be con¬ 
sidered as fully external and independent: think of the case of a simultaneous 
death caused by a car accident in consequence of the fact that the one of the two 
that was driving had an heart attack. Clearly there are factors affecting both 
the simultaneous and the separated deaths: health care, social and government 
assistance, affective relationships and so on. 

The idea of a common factor affecting all the components of a random vector 
is at the basis of the dependence structure represented by Archimedean copulas. 
Such a dependence can be in fact obtained starting from a vector of i.i.d. expo¬ 
nentially distributed components and dividing all of them by a positive random 
variable representing a common stochastic intensity (see Marshall and Olkin, 
1988). This is the case of Archimedean copulas with completely monotone 
generator. The case of /c-dimensional Archimedean copulas with fc-monotone 
generator can be similarly constructed starting from a vector uniformly dis¬ 
tributed on the unit fc-dimensional simplex and multiplying each component by 
a non-negative random variable, representing again a common factor affecting 
all the random vector (see McNeil and Neslehova, 2009). In any case, all the 
random variables involved are modified by the common factor in the same way 
and, as it is well known, the induced dependence is symmetric. 

In this paper we generalize the Li and Pellerey (2011) setting allowing for 
a dependence structure of Archimedean type among the original shocks arrival 
times. The restrictive fact of considering an exchangeable dependence struc¬ 
ture represents a first step in the perspective to include dependence: this is a 
compromise between mathematical tractability and realism of the assumptions. 

More precisely, we assume for (Ai, A 2 , A 3 ) a joint survival distribution func¬ 
tion more general than that in ([T]), that is 

Fix, y, z) = GiHi{x) + H 2 {y) -K Hs{z)) 

where G : [0, -l-oo) —>• [0,1] is the generator of a three-variate Archimedean cop¬ 
ula: in this case the associated copula is Archimedean with generator G and 
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Hi{x) = G ^(Fi(a;)) where Fi the marginal distribution of Xi. The correspond¬ 
ing survival distribution of the random vector {Mi, M 2 ) is 

F{xi,X 2 ) = G{Hi{xi) + Fl 2 {x 2 ) -b i?3(max(a:i,a:2)) 

and we call it Archimedean-based Marshall-Olkin distribution and the copula 
that it induces Archimedean-based Marshall-Olkin copula. 

As we will see, the family of the copulas so generated contains the two- 
dimensional Archimedean copulas, the Generalized Marshall-Olkin copulas and 
other well known families as specific cases. More precisely, we will show that 
any Archimedean-based Marshall-Olkin copula can be obtained through a suit¬ 
able, in general asymmetric, distortion of a bivariate Archimedean copula with 
generator G. In the case of symmetric distortions, we recover, even if under 
more restrictive assumptions, the generalization of Archimedean copulas intro¬ 
duced in Durante et al. (2007). In the case of linear distortions, we recover a 
proper subset of the Archimax copulas of Caperaa et al. (2000) that, when G is 
the Laplace transform of a positive random variable, coincide with the bivari¬ 
ate Scale Mixture of Marshall-Olkin copulas studied and applied in Li (2009), 
Bernhart et al. (2013) and Mai et al. (2013). 

The impact of the distortion on the dependence structure induced by an 
Archimedean-based Marshall-Olkin copula is analyzed by studying its Kendall’s 
function, its Kendall’s tan and its tail dependence parameters. 

The paper is organized as follows. In Section[5]we introduce the Archimedean- 
based Marshall-Olkin distribution. In Section [3] we derive the Archimedean- 
based Marshall-Olkin copula. Section |4] is devoted to the analysis of some de¬ 
pendence properties of the copulas introduced: in particular, the expression 
of the Kendall’s function and of the Kendall’s tau are derived and, in some 
particular cases, the tail dependence parameters are calculated. 

2 The Archimedean-based Marshall-Olkin Dis¬ 
tribution 

Let (n, P) be a probability space and {Xi,X 2 ,X^) be a random vector such 
that P(Xi > 0) = 1 for every i = 1,2,3. We assume that their joint survival 
distribution is of type 

F{x,y,z) = G{Hi{x) + H 2 {y) F H^iz)) 


where 

• G : [0,-boo) —>■ [0,1], with G(0) = 1 and G' exists on (0,-boo), it is non¬ 
positive, non-decreasing and concave; if xq = inf{a; > 0 : G{x) = 0}, we 
have that G is strictly decreasing on [0,a;G) and so its inverse function 
G~^ : (0,1] —>■ [0,a:G) is well defined: we can extend G~^ to the whole 
interval [0,1] by setting G“^(0) = xq- 
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• for every i = 1,2,3, Hi is continuous and strictly increasing for a: > 0, 
Hi(x) = 0 for a; < 0 and lim Hi{x) = xq- 

X^-\-OQ 

Under these assumptions, the joint survival distribution can be rewritten as 

nx,V,z) = G{G-\F^{x)) + G-\Uy)) + G-\Uz))) ( 2 ) 

with Fi = G o Hi. Thanks to Sklar’s theorem and Theorem 2 in McNeil and 
Neslehova (2009), the Fi, i = 1,2,3, are the marginal survival distributions of 
the random variables Xi and the dependence structure of the random vector 
{Xi,X 2 ,X 3 ) is of Archimedean type with generator G. 

Clearly, this setting includes, as a specific case {G{x) = e“®), the family 
of the Generalized Marshall-Olkin distributions introduced in Li and Pellerey 
( 2011 ). 

In what follows, with abuse of notation, we will denote with /(+oo) the 
lim f{x), when it exists. 

ai—>-+oo 

Remark 2.1. In Muliere and Scarsini (1987), the particular case with G(x) = 
e~^ and Hi{x) = XiH{x), Xi > 0, is studied and justified. 

The choice of functions of type Hi(x) = XiH(x) is allowed for anyG provided 
that XG = + 00 . If XG < + 00 , the requirement Hi{+co) = xg is satisfied if 
and only if Xi = X 2 = A 3 which corresponds to the case in which the vector 
{Xi, X 2 , X 3 ) is exchangeable. 


Let us now consider the random vector (Mi, M 2 ) defined by Mi = min(Ai, A 3 ) 
and M 2 = min(A 2 ,A 3 ). The corresponding survival distribution function, for 
ti,t 2 > 0 , is 

FMi,M2itl:t2) = P(Mi > ti,M2 > ^ 2 ) = 

= P(Ai > ti, A2 > t2,^^3 > max(ti,t2)) = ( 3 ) 

= G{Hi{ti) + H2{t2) + i^3(^lax(^l,^2)))■ 

We call the above distribution Archimedean-based Marshall-Olkin distribution. 
Setting Ki(ti) = Hi{ti) + H^lfi), the marginal survival distributions are 


FM,{ti)=G{K,{U)),i = 1 , 2 . (4) 

Remark 2 . 2 . Notice that if xg < + 00 , then .Fmi,M2(U, ^2) = 0 on {(ti,t2) : 
iLi(ti)+472(^2) + M3 (max(ti,t2)) > xg} and FuAti) = 0 for all U > K~\xg)- 
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2.1 The singular component 

Considering insurance contracts based on the lifetimes Mi and M 2 , it is of some 
relevance to separate the impact on the joint distribution (and, equivalently, 
on the price) of simultaneous failure or default or death (depending on the 
application in analysis) from the separated ones. 

The simultaneous failure is clearly identified by the singular component of 
the distribution. It is in fact well known that Marshall-Olkin distributions 
(as well as the Generalized Marshall-Olkin distributions introduced in Li and 
Pellerey, 2011) admit a singularity: clearly this fact continues to hold in our 
extended setting. 


To simplify the notation, we set 
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iL(t)=^iL.(t). 

i=l 

If M = min(Xi,X 2 ,X 3 ), its survival distribution function is given, for t > 0, 
by 

P(M >t) = ¥{Xi> t, X 2 >t,X3>t) = G{H{t)). 

Proposition 2.1. Assume that G is twice differentiable and that each Hi is 
differentiable on (0,-|-oo). Then 


where 

and 


Fm.mM = <t,M2<t) = 

= 1 + G{H{t)) - G{Ki{t)) - G{K2{t)) 

= - f H’^{x)G'(H{x))dx. 

Jo 


( 5 ) 


Proof. ([S]) is trivial. 
Since 
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-pAMo(h,0) = 


G" {Hi{ti) + Lf2(0)) F['i(ti)K2(t2) if ti < t2 
G" {Ki{ti) + iL2(0)) K'i{ti)H!2{t2) if h > O 


dt2dti 

it is a straightforward computation to show that 

P(Mi <t,M2< fM2 > Ml) = G{H{t)) - G{K2{t)) - H'i{x)G' (i?(a;)) dt 
and 

P(Mi <t,M 2 < t,M 2 < Ml) = 1 - G{Ki{t)) + K[{x)G' (m(x)) dt 
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and 


<t,M 2 < t,Mi = M 2 ) 



h;,{x)g' (i?(x)) 


dt. 


□ 


As a consequence 


P(Mi = M2) = - 


r+oo 


/O 


f-H ^(xg) tt! 


dt 

H;i{t)G' dt 

it)G' [H{t)) 


= E 


mt) frr 

Tm” 

M'(M)" 


dt 


H'{M) 

Example 2.1. 1. If G{x) = e~^, then 


and 


= J exp (--ff(O) dt 


P(Mi = M 2 ) = f H'^{t) exp {-H{t)^ dt 


and this is the case of the Generalized Marshall-Olkin distributions (see 
Li and Pellerey, 2011). 

2. If XQ = +00 and Hi = XtH, then, if X = X]i=i 

F!4,,M2it) = T^^-G{XHit)) 

A 

= yP(M < t) 


and 

P(Mi = M 2 ) = ^. 

A 

independently of G. In particular, this case includes, when G(x) = e~^, 
the Marshal-Olkin type distribution introduced by Muliere and Scarsini 
(1987), and, if Hi(x) = XiX, the standard Marshall-Olkin case distribu¬ 
tion. 
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3. More in general, let us assume that H 3 is proportional to Hi + H 2 , that is 
H 3 = c{Hi + H 2 ), which is equivalnet to: -Finin(Mi,M 2 )(i) = G [\H 3 {t)). 
If xa < +c», necessarily c = ^, while any choice of c > 0 is allowed if 
XG = + 00 . In any case, we have 

= -^P(M < t) 

and 

P(Mi = M2) = 

c+l 

that doesn’t depend on G. 

4 . Both the facts that is proportional to the cumulative distribution 

function of M and that the singular component is independent of G con¬ 
stitute a peculiarity of the above case. 

IfG{x) = {x-\-l)~^/^ (Clayton generator) with 9 > 0 and Hi{x)-\-H 2 {x) = 
H^{x) + H 3 {x) then 

(1 ~ (1 + ^3(0) “ 

but 

P(M < t) = I - {1 + Hsit))-^. 

Moreover,this is equivalent to 

P(Mi = M2) = ^ 

that depends on 6 . 

Similarly, if Hi{x) + H 2 {x) = e^^G) _ H^ix) — 1 then, again, 

with 

P(M <0 = 1- 

and 

P(Mi=M2) = ^. 

3 The Archimedean-based Marshall-Olkin Cop¬ 
ula Function 

If Ui = G{Ki(ti)), from (|4]), we get ti = K~^{G~^{ui)). Substituting it in 
m, thanks to Sklar’s theorem, we obtain the copula associated to the vector 
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{Ml, M 2 ), that is 


Cmi,M2{ui,U2) 


Setting Di{x) 


^ ^ r GiHioK^\G-^iui)) + G-^{u2))iiK-\G-\ui))<K-\G-^{u2)) 

\ G{G-Hui) + H2oK^\G-Hu2))) a K^\G-^{ui)) > K^\G-\u2)) 

( 6 ) 

Hi o ^{x) for X £ [0, xg] and z = 1, 2 we have 


C'mi,M2(ui, U 2 ) 


G {Di{G-\ui)) + G-\u2)) if Ki\G-\ui)) < K^\G-\u 2))) 
G {G-\ui) + D2{G-\u2))) if K^\G-\ui)) > K^\G-\u2))) 


Let us analyze the set 

F = {{ui,U2) £ [0,1]2 : K-\G-\ui)) = K-\G-\u 2 ))} = 

= {{ui,U2) £ [0, If : Di{G-\ui)) + G-\u2) = G-\ui) + D2iG-\u2))}. 

(7) 

We have 


• if Di{xg) > D 2 {xg), then [0,iLf ^(xg)] ^ [0,iL^^(a;G)] and the set ([T]) 

is the graph of the function U 2 = h{ui) = G o K 2 o o G~^{ui) for 

Ml G [0,1], where h is strictly increasing, /i(l) = 1 and h{0) = P{M 2 > 
K^\xg))>0-, 

• if Di{xg) < D2{xg), then [0, iLf ^(xg)] ^ [0, iL^^(a;G)] and the set ([7]) 

is the graph of the function mi = h{u 2 ) = G o Ki o o G“^(m 2 ) for 

M 2 G [0,1], where h is strictly increasing, h(l) = 1 and h{0) = P{Mi > 
K-\xg))>0-, 

• if Di{xg) = D 2 {xg), then the set d?]) can be represented through both h 
and h that are strictly increasing, h{l) = h{l) = 1 and h{0) = h{0) = 0. 


It follows that 


• if L>i(a:G) > D2{xg), 


C'mi,M2(wi,M2) 


G (Di{G-fui)) + G-fu 2 )) if m 2 < h{ui) 
G (G-i(ui) + D 2 {G-\u 2 ))) if m 2 > h{ui) 


( 8 ) 


• if L>i(a;G) < D2{xg), 


Gmi,M2(wi,M2) 


G (Di(G-1(mi)) + G-1(m 2)) if Ml > h(M2 ) 
G (G-1(mi) + D2{G-\u2))) if Ml < h(M2) 


(9) 


Clearly, we recover a copula of type ([S]) from a copula of type (|S]) by inverting 
Di with D 2 and viceversa. 

Remark 3.1. Notice that Dfx) and x — Dfx) are eontinuous and strietly 
increasing on [0,a;G]- Moreover 

Di{x) < X, for X > 0 (10) 

and, if XG = + 00 , lim Dfx) = +00 and lim x — Dfx) = + 00 . 

X—¥Xg x—^xq 
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In next result we will prove that for any suitable pair of functions Di and 
Z? 2 , the copula in ([8]) and can be obtained starting from a three-variate 
distribution function of type (H)) and following the construction presented above. 

In what follows we set Di{x) = x — Di{x). The following result holds: 

Lemma 3.1. Let Di^D^ be two strictly increasing and continuous functions 
on [OjXg] with Di{0) = 0, Di{x) strictly increasing and such that, if xq = 
- 1 - 00 , lim Di{x) = Di{xc) = +oo and lim Di{x) = Di{xG) = -l-oo. Then, 

X—^XQ X—¥Xq 

for every defined on [0,-|-oo), continuous, strictly increasing and such that, 
HsiO) — 0 and H^^+oo) = xq, the function Hi(x) = D~^{H^ix)) — H 3 {x), for 
X G [0, is the unique solutions of 

Di{x) = HiO {Hi + H3)~^{x) for all x G [0,a:G]. (11) 

Proof. Notice that Hi{x) is strictly increasing on [0, Hf^{Di{xG))] and satisfies 
Hi{ 0 ) = 0, Hi{Hf^{Di{xG))) = Di{xG) and it is the only solution of 


Di{Hi{z) + H3 {z)) = H3{z) 

z€[0,H3\Di{xg))] 


which is equivalent to 

Di{Hi{z) + H3 {z)) = Hi{z) 

z G [0,iJ3'^(A(a;G))] 

(12) 


Setting X = Hi{z) + H 3 {z) G [OjXg] in ([TH) we get (fTTI) . 

If xg < - 1 - 00 , Hi can be clearly extended so to be defined on the set [0, -l-oo), to 
be strictly increasing and with Hi{+oo) = xg- 

□ 

Remark 3.2. Notice that Di and Di satisfy the same assumptions and if, given 
H 3 , Hi solves 177]] . thenHi{x) = D~^{H 3 {x))—H 3 {x) = {d~^ — Id^ {H 3 {x)) 

(where Id is the identity function) solves Di{x) = Hi o {Hi -1- H 3 )~^{x) for all 
X G [0,a:G]- 

Example 3.1. 1. If Di{x) = aiX, oti G (0,1), i = 1,2, then Hi{x) = 

T^H 3 {x), fori = 1,2 andxG [o,R 3 "^((l - ai)xG)]. 

2. IfDi{x) =x-^^ log{aiX+l), ai > 0, i = 1, 2, then Hi{x) = ^ (^e°‘iH 3 D) - l)- 
H 3 {x), fori = 1,2, and x G [0, R 3 "^(^ log(Q:ia;G + !))]■ 

3. If Di{x) = X — {\/aiX -I- 1 — 1 ), a.i > 0, i = 1,2, then Hi{x) = 

aiH 3 {x) -I- H 3 {x), for i= 1,2 and x G [0, Hf^{-^{3faiXG -|- 1 - 1))]. 

4. IfDi{x) = ^ {3/a~xTT- 1), ai > 0,i = 1 , 2 , then Hi {x) = (^i/l + iaiH 3 {x) - 1 

for i = l,2 and x G [0,Hf^{xG - ^{^UiXG + 1- 1))]. 
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Definition 3.1. Let G : [ 0 ,+c») —>• [0,1], with G'(O) = 1 and such that G' 
exists on ( 0 ,+c»), it is non-positive, non-decreasing and concave. If xq = 
inf {a; > 0 : G{x) = 0}, let Di,D 2 he two strictly increasing and continuous 
functions on [OjXg] with Di{ 0 ) = 0, Di{x) strictly increasing and such that if 
xg = + 00 , lim Di(x) = Di{xG) = +oo and lim Di(x) = Di^xc) = +oo. If 

x—>^xc x—^XG 

we set f{x) = o Di{x) (when this is well defined), we call Archimedean- 
based Marshall-Olkin copula with generator G and distortions Di , D 2 the copula 
function Camo so defined: 


• if Di{xg) > D2{xg), 


Camo{u, v) 


G{Di{G-^{u)) + G-\v)) ifv<h{u) 
G{G-\u) + D 2 {,G-\v))) ifv>h{u) 


where h{u) = G o f o G ^{u); 


• if Di{xg) < D2{xg), 


Gamo{u, v) 


G{Di{G-^{u)) + G-\v)) ifu>h{v) 
G (G-Hu) + D2iG-\v))) ifui<h{u2) 


(13) 


(14) 


where h{v) = G o f ^oG ^(r;). 


It is evident from the definition that the Camo copula is a distortion of the 
Archimedean copula with generator G through the functions Di. More precisely, 
if Di 7 ^ D 2 (that is if Hi ^ H 2 ), the Archimedean copula is differently modified 
above and below the curve F given in 0 and the obtained copula is obviously 
asymmetric. 

Conversely, \i Di = D 2 = D (that is if Hi = H 2 ) we get 


Gamo{u, v) 


G (D{G-\u)) + G-\v)) ifu<u 
gIg-'^{u) + D{G-^{v))} ilv>u 


(15) 


and the obtained copula remains exchangeable. The family defined in (1151) is 
contained in the class of copulas introduced in Durante et al. (2007), defined as 

A v) + V'(u V ri)) (16) 


with (j) : [ 0 , 1 ] —>■ [ 0 , + 00 ], continuous, convex and strictly decreasing, : [ 0 , 1 ] —>■ 
[ 0 ,+oo], continuous, decreasing and such that ' 0 ( 1 ) = 0 and 0 — 0 increasing 
in [ 0 , 1 ]: (fT5ll is trivially recovered from (IT 6 )) when 0 ( 1 ) = G“^(l) = 0 and 
0(0 = D(G-0t)). 

Example 3.2. If G{x) = e~^ 


Camo{u, v) 


ve if V < exp{—f In u)) 

yg-£> 2 (-in«) jy^; > exp(—/(—Inu)) 


and we recover the Generalized Marshall-Olkin copula introduced in Li and 
Pellerey (2011). In particular, if, for i = 1,2, Di{x) = aiX with Ui G (0,1), we 
recover the classical Marshal-Olkin copula. 
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Example 3.3. If Di{x) = UiX for i = 1,2 and ai G (0,1), then 


Camo{u, v) 


G (aiG-i(«) + G-i(z;)) ^f y < G (u)') 

G {G-\u) + a^G-Hv)) ifv>G \}^G-\u)) 


(17) 


This particular specification constitutes a subclass of the family of copulas called 
Archimax copulas introduced in Caperad et al. (2000). Archimax copulas are 
defined as 


GgAu,v) = G ({G-\u) + G-\v)) A 


G-\u) 


G-i(w) + G-i(z;) 


(18) 


where A : [0,1] —>■ [ 1 / 2 , 1 ] is a convex function such that max(l, 1 — t) < Aft) < 1 
for all t G [ 0 , 1 ] and copulas in can be obtained choosing 


Ait) 


iai-l)t + l ift < 
(1 - a2)t + a2 ift> 


The same fact holds if one considers piecewise linear distortions Di, to which 
correspond piecewise affine functions A in representation m- 

Moreover, the subset of copulas of type when G is the Laplace trans¬ 
form of a positive random variable, coincides with the family of bivariate Scale 
Mixture of Marshall-Olkin copulas (see Li (2009), Bernhart et al. (2013) and 
Mai et al. (2013)). These copulas are obtained as the copulas associated to the 
random vector (■§)",-f^) where {Zi,Z 2 ) has a bivariate Marshall-Olkin distri¬ 
bution and Y is an independent positive random variable. In the exchangeable 
case. May et al. (2013) propose an alternative construction: let £ 1,62 be i.i.d. 
unit exponentially distributed random variables, M > 0 be a random variable 
with Laplace transform G{x) and At 7 ^ 0 6 e a Levy subordinator with Laplace 
exponent 4'; assuming that they are mutually independent, define 


Tk = infji > 0 : Mt > fife} 

where Mt = Atv 7 G-i(i-p(t))/«'(i) vi^) ® given distribution function. The 
copula associated to {ti,T 2 ) is 

G{u,v) = G I^G~^(min(u,v)) + G“^(max(u, 

which is obviously of type |17| ) with ai = a 2 = ’ 


Under the assumptions of Proposition 12.11 


(<H ^( xg ) 


iJ^(t)G' (i7(t)) dt = F{Mi 


M 2 ) = 


= PiFf(,l{U) = Ff(,l{V)) = 

= FiKf^ o G{U) = o G(U)) = 
= F{F). 
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Hence, if G is twice differentiable and each Di is differentiable on [0 ,xg], the 
copula function Camo{u,v) has a singularity on the curve F (see ([7])), whose 
mass, expressed in terms of the distortions, is 

rT-\xG) 

¥{F) = - G'{T{x))dx 

Jo 

where 

T{x) = Di^{x) + D 2 ^ix) — X (19) 


It is a known fact that two Archimedean copulas with generators Ga and 
Gb coincide if and only if there exists a > 0 such that Ga{x) = Gb{o-x) (see 
Corollary 2.2.6 in Alsina et ah, 2006). In what follows we will present the 
analogous result for the Archimedean-based Marshall-Olkin copula functions. 

Theorem 3.2. Let G^ and G^ he two Archimedean-hased Marshall-Olkin cop¬ 
ula functions with generators Ga o,nd Gb and distortions Ai, A 2 and Bi, B 2 , 
respectively. 

CA ^ fjB 

if and only if there exists m > 0 such that 


Gb{z) = GA{inz) for z G [0, + 00 ) 


( 20 ) 


and 

Bi{z) = —A^{mz) for z G [0,xgb] , * = 1,2. (21) 

m 

Proof. Since we are looking for G^ = G^, necessarily the corresponding singu¬ 
larity sets F in (EU must coincide. Hence and must be of the same 
type ra or m- We show the proof in case (ITSl) being the other perfectly 
identical. Hence, we assume Ai{xga) ^ ^ 2 ( 3 ^ 0 ^), Bi{xgb) ^ B 2 {xgb) and 


G^{u,v) = 


Ga (Ai(G^^(u)) -b G^^(w)) if u < hA{u) 
Ga {G^\u) + A 2 {G^Hv))) ■Av>hAiu) 


and 

gB(b b) = 1 {B,{Gi\u)) + Gl\v)) if < hB{u) 

[U,V) \ Gb {g-b\u) + B2{G-s\v))) \iv>hB{u) ■ 

Let h{u) = min(/i/i(u),/ ib(u)). For G^ = G^, necessarily 

Ga {Ai{G-/{u)) + G)(1(z;)) = Gb (Hi(Gs'(w)) + Gs'(^)) (22) 


for u G [0,1] and v < h(u). We set if = G^^ o Gb '■ [0,xgb] [0,xg.4], 

X = Gg^{u) G [0,xgb] and y = Gg^{v) G [0,xgb]- So (1^ can be rewritten as 


Ai {^f{x)) +'ij}{y) = ifi^Bi (x) + y) (23) 
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with y G [OjXGb] and x < G^^(/i“^(Gb( 2/))). Since Ga and Gb are differen¬ 
tiable on (OjXG^) and (OjXGb)) respectively, and their derivatives are strictly 
negative, ijj is differentiable on (0 ,xgb) as well. After differentiating (1^51) with 
respect to y, we get 

i’'iy) = ’’P'iBiix)+y). 

This holds whenever 0 < x < Gg^{h~^{GB{y)))', hence being locally con¬ 
stant, is also globally contant: then there exists m such that ip'{y) = m for all 
y G (0,xgb)- It follows that ipiv) = for y G [0,xgb]i and since ip is 

strictly increasing, necessarily, m > 0. Moreover, since ip{0) = 0, q = 0. 

Hence 

Gb(x) = GA^mx), X G [0,xgb] (24) 

and xgb = 

Substituting this relation in (1^ we get 

Ai(G^i(n))=mi?i 

and, setting z = S [O, -^], we obtain 

Bi{z) = —Ai(mz). 
m 

In order to guarantee that 

Ga {G-j,\u) + A2{G-/{v))) = Gb (Gs'(n) + B^i^Gs^v))) 

must hold for v G [0,1] and u < mm(h^^(v), (v)), as well. Using (l24l) we 

obtain 

A,(G^\v))=mB, 

and, setting z = G [O, , we obtain 

B2{z) = —A2(mz). 


The converse is trivial. 


□ 


4 Dependence Properties 

In this Section we will consider the concordance measure induced by the Gamo 
copula by calculating the corresponding Kendall’s function and Kendall’s tau. 
Since, as we noticed, in general, the Gamo copula is an asymmetric distortion 
of an Archimedean copula, we investigate the effect of such a distortion on the 
tail dependence parameters. 

Let us start by observing that, thanks to (US, 

Gamo{u,v)>G{G-\u) + G-\v)), (u,z;) G [0,1]2. (25) 
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4.1 The Kendall’s function 

We remind that the Kendall’s function of a copula C, is the cumulative distri¬ 
bution function of the random variable C{U, V) with respect to the probability 
induced by C (see Nelsen, 2006), that is the C-measure of the set 

At = {(m,u) e [0,1]^ : Camo{u,v) < t}. 

Let Kelt) be the Kendall’s function of the Archimedean copula with generator 
G and Kamo be the Kendall’s function of the copula Camo- Clearly, from (El]) 

KAMo{t) ^ Kcit), te[0,l]. 


More precisely. 

Theorem 4.1. 


KAMoit) = Kait) + G'{G-\t)) ■ T-\G-\t)) 
where T is defined in m- 

Proof. We follow the same ideas and spirit of the proof of the analogous result 
in the Archimedean case (see, for example, Theorem 4.3.4 in Nelsen, 2006). 

We will prove the result in the case Di{xc) > D 2 {xg), being the alternative 
case perfectly analogous. 

We start with looking for the intersection of the t-level curve {(u,u) € [0,1]^ : 
Gamo{u,v) = t}, for t > 0, and the graph of the function h of Definition 13.11 
that is we solve for u the equation 

G(Di(G-i(u)) + G-\h{u))) = t 

from which 

Di(G-i(u)) + /(G-i(u))=G-i(t) 
or, setting g{x) = Di{x) + f{x), 

giG-Hu)) = G-\t). 

Since, thanks to the assumptions, g is invertible, the unique solution is 

Ut=Gog-^ {G~^{t)) €[t,l] (26) 

and, if 

vt=h{ut) =Go f o g-^ [G~^{t)) , (27) 

(ut, Vt) is the unique intersection of the level curve and the graph of h. 

We split the set At into three regions: the rectangle Rt = [0, Ut] x [0, Vt] and 
the sets 


Bt = {{u,v) G [0,1]2 : u G [ut, l],u < G{G-\t) - D,{G-\u)))} 
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and 


B2={{u,v) G [0,1? ■■ve[vt,l],u<G{G-^{t)-D 2 {G-\v)))} 


Obviously 

P(i?t) = G(ut,vt) = t. 

Let us compute now P(i?i) and P(i? 2 )- We start with P(-Bi). 

Consider the partition of the interval [ut, 1] given by the points 

tk = Go {G-?ut)) (1 - . 

Notice that to = ut < ti <■■■< tn = 1 and that, since G o is uniformly 
continuous on [0,£>i(G'“^('Ut))], the width of the partition so defined goes to 
zero as n —>■ +oo. 

For M G [«-(,!], the t-level curve in Bi is 

viu) = G {G-?t) - Di{G-\u))) . 

For fc = 1,... ,n, let Rk = [4-1,4] x [0, v(4-i)]- Clearly 

P{Rk) = G{tk,v(tk-i)) - G{tk-i,v{tk-i)) = G(tk,v(tk-i)) - t 


and 

G{tk,v{tk-i)) 


If Sr, = ELi Rk, 


G {D^{G-?tk)) + G-\v{tk-i))) = 

G(G-\t)-^l^^^:^\. 


P(^„) = n 


G G-^{t) 


D^{G-\ut)) \ - 

n ) 


G (G-i(t) - 

Di(G-Hut)) 

n 


Di{G-?ut)) 


and 

P(4i)= hm P(,S„) = -G'(G-i(t))i?i(G-^K))- 

n—>-+oo 


For the set B 2 , considering the partition of the interval [vt, 1] given by the 
points 

Sk = G {d 2 {G-?vt)) (1 - ^))) . 

exactly as done for Bi, we get 

P(42) = -G'{G-\t))D2{G-\vt)). 
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It follows that, using ([26|) and (l27ll . 


F{At) = P(i?t) + P(Bi) + P(B2) = 

= t-G'{G-\t)) [D,{G-\ut)) + D2{G-^vt))] = 

= t - G'{G-\t)) [D,{g-\G-\t))) + D2 o f {g-\G-\t)))\ = 

= KG{t) + G\G-\t)) [G-i(t) - D^{g-\G-\t))) - o f {g-\G-\t)))\ 

= Koit) + G\G-\t)) [Id -D^og-^-D^ofo g-i] o G-\t) 

where Id is the identity function and Kait) = t — G'{G~^{t))G~^(t) is the 
Kendall’s function of the Archimedean copula with generator G. 

But, since g{x) = T o -Di(x), where T is defined in (fTOl) . we have 

Id - DiO g-^ - D20 f o g-^ = {g - Di- D20 f)o g-^ = 

= if - D20 f)o g-^ = 

= £>20/05-1 = 

= £>105-1 = 

= T-i 

and so 

KAMoit) = Kcit) + G'(G-i(t))T-i(G-i(t)). 

□ 

In terms of the generating functions Hi, II2, II3 the Kendall’s function can 
be rewritten as 

KAMoit) = Kcit) + G'(G-I(t)) • H3 O i/-i(G-i(t)). 


In particular, if H 3 = c{Hi + H 2 ), with c > 0 (see Example 12.11) we get 

KAMoit) = Kcit) + ^G'(G-i(t))G-i(t) (28) 

c + 1 

and this case includes the case of linear distortions of Example 13.11 

Notice that (|28)) can be rewritten as KAMoit) = t — ■^G'iG~^{t))G~^it), 
and if we consider two different generators Ga and Gb so that KcAit) < 
Kcsit), but identical proportional parameter c, then KcAMoit) ^ KcgMoit)^ 
where with Kq^mo we denote the Kendall’s function of the Archimedean-based 
Marshall-Olkin copula with generator Gj, for j — A,B. 

This fact doesn’t continue to hold in general as shown in next example. 
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Example 4.1. Let Di{x) = D 2 {x) = 1 + x — y/1 + 2x. It follows that 


KAMoit) = t - ^G\G-\t)) (v'l + 4-G-i(i)- l)" 

If we consider the Frank generator Gf{x) = —^ In (l + — 1)) with pa¬ 

rameter 9 = 4 and the Gumbel generator Gg(x) = exp (— with parameter 
') = 2we have 

KGpiO.3) = 0.497 > 0.480 = ^ 00 ( 0 . 8 ) 

while 

KgfMo{0-^) = 0.341 < 0.380 = KGaMoiO-^)- 

Example 4.2. Let us consider Di{x) = x — -^{\/aiX + 1 — 1) (see Example 
\3.1\] . Then 

if.Mo(t)= KAt )+. (\/i+i(o.+«)G-.(o-1) 


4.2 Kendall’s tan 

It is known that the Kendall’s tau of a bivariate copula G is given by 

Tc = 4E[G{U, I/)] - 1 = 3 - 4 [ K{t)dt. 

Jo 

Hence, we get 


tamo = 3 — 4 / KAMo{i)dt = 

Jo 

= 3 - 4 / (Kcd) + G'{G-\t)) ■ T-\G-\t)))dt = 

Jo 

= tg-4 [ G'{G-\t)) ■ T-\G-\t))dt = 

Jo 

= tg+4 / {G'{x)) ■T~^{x)dx 

Jo 

where tg is the Kendall’s tau of the Archimedean copula with generator G. 
Obviously, as expected from (|25l) . 


tamo > tg- 


In terms of the generating functions ili, 772, Ho the Kendall’s tau can be rewrit¬ 
ten as 


tamo =tg+4 / (G'(a:)) • H^, o II~^{x)dx. 


Jo 
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Example 4.3. IfG{x) = e ^ we recover the case studied in Li and Pellerey(2011) 
and we get 

p-\-oo 


Tamo =4: / e ^{x)dx. 


Jo 


Example 4.4. Let us consider the Clayton case, that is G{x) = {x + 1) 
with 9 > 0 . 


1. If H 3 = c{Hi + H 2 ), 0 0, then 

tamo =rG + ^^ = rG + nMO^ 


where tjmo = the Kendall’s tau when G = e Notice that this 

result continues to hold for 6 G (—^, 0 ) and c = i. 

2. If L)f^{x) + Dif^ix) = e~ — 1 -\- X, with 1 > 7 > 0 (this corresponds to 
IIi{z) + Il 2 {z) = e T —Il 3 {z) — 1), we get (^Df^{x) + Df^{x) — = 

7 ln(x + l). Then 


tamo = TO 


47 

(2 + 0)2 


9^ + 29 + 47 
(2 + 0)2 


3. If Di ^(x) + £>2 ^(x) = — 1 + X with a>l, 1>7>0 

(this corresponds to IIi{z) + Il 2 {z) = + 1^ — 1 — H 3 {z)), we get 

(^Df^lx) + L)f^{x) — x'j = 7 ^(x + l )“—Then 


Tamo = tg 


47 

{2 + 9){a{2 + 9) - 9) 


1 

9 + 2 



47 


a{2 + 9)-9 


4 . If Di{x) = x—^{\/aiX + 1—1) (see Example \4-2{ ), (jj^ ^(x) + £>2 ^(x) — x^ 

2 (ai+a 2 ) ((y^l + |(ai + Oi 2 )x - . If cti + a 2 = j, we recover the previ¬ 

ous case with 7 = § and a = 2. Hence 

36»2 + 120 + 8 

~ 3(2 + 0)(4 + 0)' 

Exactly as for the Kendall’s function, the presence of distortion functions 
strongly influences the Kendall’s tau. In Table[I] the values of the Kendall’s tau 
of some Archimedean copulas and the corresponding values of the Archimedean- 
based Mashall-Olkin copulas with same generator but distortions Di{x) = 
D 2 {x) = 1 + X — \/l + 2x are reported in order to illustrate how the pres¬ 
ence of the distortions can induce an inversion in the order of the concordance 
measure. 
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Table 1: Comparison of Kendall’s tau value with and without distortions. 


Generator 

Parameter 

XG 

XAMO 

Clayton 

9 = 2 

0.5 

0.78539 

Gumbel 

0 = 1.8 

0.375 

0.81636 

Frank 

II 

0.388 

0.906 


4.3 Tail dependence 

We recall that the upper and lower tail dependence parameters of a copula C 
are given, respectively, by 


Xjj = 2 — lim 
nfl 


1 — C{u, u) 
I — u 


Cit 


and Xl = hm- 

uiO u 


when these limits exist and that if Xu G (0,1] (A^ G (0,1]) we have that C has 
upper(lower)-tail dependence (we refer to Section 5.4 in Nelsen (2006), for more 
details). 

Let Xf^^^ and Xu^^ be the lower and upper tail dependence parameters 
of the copula Camo and and A^ be the tail dependence parameters of the 
corresponding bivariate Archimedean copula with generator G. 

Let Ai = {x G [0 ,xg) : Di{x) > D 2 {x)} and A 2 — A^. Notice that 
G~^(u) G Ai if and only if {u, u) lies on or below the curve F defined in ©• 


4.3.1 Af^*^ 

If XG < + 00 , obviously Xf^^^ = 0. 

Let us now consider the case xq = + 00 . We have 


Camo{u,u) = G {G-\u) + Di{G-\u))1aAG-\u)) + 7A2(G-1(m))1^,(G-1(u))) . 
Setting X = G~^{u), 


^AMO _ jj^C'amo(u, u) 


u 4 ,o u 

G{x + Di (x) 1^1 (x) + B 2 (x) 1 a2 (x)) 


= lim 

X—>■+00 


G{x) 


If there exists x > 0 such that {x, + 00 ) C Ai then 

XfMO = lin, G(a: + A(x)) ^ 
rc-f+oo G(x) 

The value of the above limit clearly depends on the type of decay of G to 
zero. If the decay is of polinomial type, that is there exist c, 7 > 0 such that 
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G{x) ~ cx ^, then 

a:—)-+oo 


lim 

X—>- + oo 


G(x + Di{x)) 

W) 


= lim 1 + 


Di{x) 


Hence, if lim = fl g [Q 1] exists, we have that 

x^+oo “ 

^AMO ^ (l+;3)-7_ 

If the decay is exponential of type G{x) ^ ce~°'^^, with a,c,^ > 0, we have 

21—>- + 00 

G(x +A(x)) ^ ^_aax+D{x)r-x-’) ^ 


lim 

a;^+oo G(X 


a:—>-+cxD 


= lim e—^((1+"^)^-^). 


a:—>-+oo 


If lim =/3 g (0,1], then =0. 

ai—>-+oo ^ 

If lim = 0, then 

ai—>-+oo ^ 


ai—>-+oo ai—)-+oo 

and, if 7 > 1, then X^^^ = 0 while, if 7 < I, then 

Di(,x) 

XAMO ^ 

Example 4.5. Assume that there exists a; > 0 such that (x, + 00 ) C Ai with 
lim = Pi. lEe have 

ai—>-+oo ^ 

• in the Clayton case G{x) = {x + 1)“^/®, 9 > 0, = (I + Pi) > 

2-V« = Ag; 

• in the Gumhel case G(x) = e~^ ^ , 9 > 1, with pi G (0,1], = 0 = 

Ag; 

• in the Frank case, G{x) = In (l + — 1)), 0 G M, X^^^ = 0 = 

Af. 


4.3.2 X§^^ 

Similarly as done for Xf^^^ we can calculate X^^^ in some specific cases. 
If there exists a; > 0 such that (0, x) C Ali then 


\AMO 


= 2 - 


1-G(x + A(a;)) 

hm-—— - 

a:-i -0 I — G{x) 
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If 


for some c, 7 > 0 , then 


1 — G(x) ~ 

xiQ 


. l-G{x + Di{x)) 
" 1 ) 1 - G(x) 


lim 


= lim I 1 + 




Di{x) 


Hence, if lim !Me 1 = p. g [0,1], then 


XAmo ^2-{1 + Piy. 


Example 4.6. Assume that there exists x > 0 sueh that (0,x) C Ai with 
lim £ [Q^ X]. We have 

a:^0+ ® 

• in the Clayton case G{x) = (x + 1)“^/^, 0 > 0, = 1 — /3i > 0 = ; 

• in the Gumbel case G(x) = e~^ ^ 0 > 1, = 2 — (1 + > 

2 - 2^/® = Ag; 

• in the Frank case, G{x) = —| In (l + — 1)), 9^0, = 

1-ft >0 = Ag. 
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